NASA/CR— 2003-212306 



Xenon Sputter Yield Measurements for Ion 
Thruster Materials 


John D. Williams, Michael M. Gardner, Mark L. Johnson, and Paul J. Wilbur 
Colorado State University, Fort Collins, Colorado 


May 2003 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/CR— 2003-212306 



Xenon Sputter Yield Measurements for Ion 
Thruster Materials 


John D. Williams, Michael M. Gardner, Mark L. Johnson, and Paul J. Wilbur 
Colorado State University, Fort Collins, Colorado 


Prepared for the 

28th International Electric Propulsion Conference 

cosponsored by the Centre National d'Etudes Spatlales (CNES), European Office of Aerospace 
Research and Development of the USAF, Snecma Moteurs, and European Space Agency (ESA) 
Toulouse, France, March 17-21, 2003 

Prepared under Grant NAG3-1801 


National Aeronautics and 
Space Administration 


Glenn Research Center 


May 2003 



Acknowledgments 


Partial financial support under Research Grant NAG3-1801 through the 
NASA Glenn Research Center is gratefully acknowledged for P.J. Wilbur and M.M. Gardner. 


This report contains preliminary 
findings, subject to revision as 
analysis proceeds. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http:/ /gltrs. grc.nasa.gov 


XENON SPUTTER YIELD MEASUREMENTS 
FOR ION THRUSTER MATERIALS 


John D. Williams,* Michael M. Gardner, Mark L. Johnson, and Paul Wilbur 

Colorado State University 
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ABSTRACT 

In this paper, we describe a technique that was used to measure total and differential sputter yields of 
materials important to high specific impulse ion thrusters. The heart of the technique is a quartz crystal monitor 
that is swept at constant radial distance from a small target region where a high current density xenon ion beam 
is aimed. Differential sputtering yields were generally measured over a full 180° arc in a plane that included the 
beam centerline and the normal vector to the target surface. Sputter yield results are presented for a xenon ion 
energy range from 0.5 to 10 keV and an angle of incidence range from 0° to 70° from the target surface normal 
direction for targets consisting of molybdenum, titanium, solid (Poco) graphite, and flexible graphite (grafoil). 
Total sputter yields are calculated using a simple integration procedure and comparisons are made to sputter 
yields obtained from the literature. In general, the agreement between the available data is good. As expected 
for heavy xenon ions, the differential and total sputter yields are found to be strong functions of angle of 
incidence. Significant under- and over-cosine behavior is observed at low- and high-ion energies, respectively. 
In addition, strong differences in differential yield behavior are observed between low-Z targets (C and Ti) and 
high-Z targets (Mo). Curve fits to the differential sputter yield data are provided. They should prove useful to 
analysts interested in predicting the erosion profiles of ion thruster components and determining where the 
erosion products re-deposit. 


INTRODUCTION 

As the interest in ion thrusters operated at high specific impulse builds, it is becoming increasingly 
important for analysts to have sputter yield data for a wide range of incident ion angles and energies. 1,2 Both 
total sputter yield (atoms/ion) and differential sputter yield (atoms/ion/steradian) are needed. These data are of 
interest for calculating erosion of spacecraft surfaces and thruster components that can be struck by high-energy 
ions produced in the beams of these thrusters. 3,4 Of equal interest is sputtering caused by ions with lower 
energies that are created via charge exchange reactions near ion thruster grid surfaces 1 or by potential hill 
structures located within ion thruster discharge chambers. 5,6 Sputter erosion of surfaces at cathode potential that 
face the interior of discharge chambers can also occur via bombardment from multi-charged ions produced 
within the discharge plasma. 7 These multi-charged ions can obtain energies corresponding to multiples of the 
anode to cathode potential difference, which can be well above the sputtering threshold 8 for most materials. 

One important use of differential sputter yield data is in numerical simulations that are needed to 
calculate erosion profiles of accelerator grid surfaces and subsequent contamination buildup on adjacent grid 
surfaces, thruster components, and spacecraft surfaces. Detailed knowledge of the evolution of the eroded 
accelerator grid surface would allow thruster performance to be predicted over long mission times. Furthermore, 
sputter yield data are also needed to describe back sputter rates onto thruster and vacuum chamber surfaces from 
beam dumps in facilities where life tests and acceptance tests are performed. These data are needed to 
determine if the test facility is affecting the test results. In response to these needs, we have initiated an 
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experimental and numerical program to provide sputter data corresponding to conditions unique to high specific 
impulse thrusters. 

In this paper, we first present a description of our measurement technique, which consisted of a quartz 
crystal monitor that was swept at constant radial distance from a small target region where a high current density 
xenon ion beam was aimed. This section is followed by a presentation of data that has been collected to date. In 
addition to experimental data, curve fit coefficients are listed in tabular form to simplify the use of the 
measurements in models of electric propulsion devices. 

SPUTTER YIELD MEASUREMENT TECHNIQUE 

Sputter yield testing was performed by first clamping a planar target to a water-cooled mounting plate in 
a vacuum system, evacuating, and then baking the system. No-flow pressures below lxlO" 6 Torr were achieved 
before testing was started. The basic test apparatus, which is shown schematically in Fig. 1, includes an ion 
source equipped with a well-focused, 2.5-cm diameter ion optics system. The ion beam is directed onto the 
center of the target surface, which is located 23 cm downstream of the ion source grids. This target surface, 
which is large compared to the beam size, can be rotated to change the angle of incidence of the beam ((3). For 
this study, incidence angles were varied from (3=0° to (3=70°. 


Back Half 

Front Half 

a > 0 


Quartz 

Crystal 

Monitor 


Fig. 1 Sketch of Test Apparatus. 

Sputtered material was sensed as a mass accumulation rate by the quartz crystal monitor (QCM), which 
is also shown in the Fig. 1. The QCM was mounted on an arm that rotated about the same axis as the one used 
for target rotation. The QCM sensor was swept on an arc that was 17.8 cm in radius and was oriented so its 
sensing face remained pointed at the midpoint of the ion beam impact zone as suggested in Fig. 1 . Mass 
accumulation rates were typically measured over the atom ejection angle (a) range from a=-90° to a=+90° 
relative to the target surface normal vector in 10° to 15° increments. The polar angle range from 0 to +90 is 
referred to as the front half of the hemisphere located above the target, and the range from -90 to 0 is referred to 
as the back half. Material sputtered in the forward direction is detected in the front half of the hemisphere and 
material sputtered backward (toward the incident ion beam for example) is detected in the back half of the 
hemisphere. In is noted that QCM measurements were not made at positions where the monitor would intercept 
the ion beam and where structural interference between the QCM monitor and vacuum chamber structure 
occurred. In addition, the quartz crystal in the monitor was replaced each time a new material was tested to 
ensure adequate bandwidth of the QCM was available for the next test. 


a < 0 


Ion 

Source 


Target 

Normal 



QCM data were obtained as changes in mass on the QCM crystal over prescribed intervals of time 
(several minutes) at each QCM position. Data collection was repeated as necessary to assure the mass 


NAS A/CR— 2003-2 1 2306 


2 




accumulation rate had reached a stable value before final measurements were recorded. During data collection, 
the total and partial pressures within the vacuum facility were monitored in nitrogen equivalent units using a 
residual gas analyzer (RGA). The total pressure was typically in the low 10' 5 Torr range and xenon was the 
dominant gas in the chamber at a nitrogen equivalent partial pressure that was very close to the total pressure. 
Nitrogen (and/or carbon monoxide) was the dominant impurity present during data collection at a partial 
pressure that was a few percent of the total pressure. Other impurities that were present during data collection 
were at an order of magnitude lower partial pressure and included oxygen, water vapor, and helium. When the 
ion beam source was turned on and ions began to strike a given target material, the 28-amu partial pressure 
signal would typically rise and then decay slowly over a period of hours, thereby suggesting that either nitrogen 
was being released from warming surfaces or CO was being produced from the hot filament used in the 
discharge chamber of the ion source. Oxygen and water vapor partial pressures would generally drop slightly 
when the beam was first turned on. No other impurities that could be detected by the RGA (<300 amu) were 
observed. 


One important concern during sputter yield measurements is the need for dynamically clean target 
surfaces. 9 This is because the presence of background gases adsorbed onto the target surface can mask the true 
sputter yield behavior of the material. A particularly troublesome background gas that is known to affect sputter 
yields of metal targets is nitrogen. One can estimate the fraction of the target covered with background gases 
using the following equation 


fin 


Pn K 


< 0.1 


( 1 ) 


In Eq. (1), p n represents the sticking coefficient of a neutral gas atom (n), n n the rate of arrival of neutral 
particles to the target surface, Y n the sputter yield of an adsorbed atom, and the rate of arrival of bombarding 

ions. As a general rule of thumb, a value of r| n less than 0.1 will ensure that background gases will not affect 
sputter yield measurements significantly, and, during our experiments, worst case estimates of the background 
gas shielding effect (r| n ) calculated using Eq. (1) were less than 0.05. 


Measurements made at various polar angles (a) using the QCM provide mass accumulation rates [R(a)] 
in gm/sec. These can be converted into units of atoms/sec by dividing the mass rates by the molecular mass of 
the sputtered material (M in gm/mol) and multiplying by Avagadro’s number (N A in atoms/mol). These results 
are converted into fractional sputter yields by dividing by the ion arrival rate, i.e. by the ion current (J) in 
coul/sec over the ion charge (q) in coul/ion. However, these yields represent only those atoms sputtered onto the 
QCM sensor through the solid angle subtended by the sensor from the region of beam impact. Because the 
sensor is swept through a circular arc and its surface always remains normal to the radius vector this solid angle 
is the same for each measurement. Assuming the ion impact region is small compared to the arc radius, the solid 
angle is equal to the sensor area (A s ) divided by the square of the arc radius (r). Dividing the fractional sputter 
yields by the solid angle subtended by the sensor gives angle-dependent differential sputter yields [y(a)]. The 
equation for differential sputter yield in atoms/ion/steradian developed as it has been described above is: 


y (a) = 


R(a) N^qr" 
M J A c 


( 2 ) 


Typical data obtained for 1 keV xenon ions sputtering titanium at a beam current of 15 mA are shown in Fig. 2. 
The constants needed in Eq. (1) to convert raw rate data [R(a)] into differential yield data for this figure are: 


N a = 6.023xl0 23 atoms/mol q = 1.6xl0' 19 coul/ion r = 17.8 cm 

M = 47.9 gm/mol (titanium) J = 15 xlO' 3 coul/sec A s = 0.535 cm 2 . 


It should be possible to obtain the total sputter yield by integrating a function describing the variation in 
differential yield with a over the full solid angle range through which atoms can be sputtered ( 2n steradians). 
For the case of normal incidence, where an axi-symmetric distribution is assumed to exist, the differential solid 
angle is given by 
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d i// = 


2nx sin (a) r da 


= 2 n sin (a) da 


The total yield is then given by 

71 71 

2 2 

Y = | 2 n sin(a) y(a) da « J n sin(|a|) y(|a|) da . (4) 

0 71 

_ 2 

For the normal incidence data in this report, the total yield was determined by fitting the differential yield data 
like those for normal incidence in Fig. 2a using up to a sixth-order polynomial and performing the integration 
analytically. 10 

In order to obtain total yield data for non-normal incidence data, it would be necessary to measure 
differential yields over the entire hemisphere above the target [i.e. obtain y(a, 4>) where 4> is the azimuthal angle] 
and then integrate the equation: 


| J 2 2 71 sin(a) y(a, §) da d 4 > 


However, approximate values of the total yields can also be obtained for non-normal incidence yield profiles 
(like the one for 30° shown in Fig. 2b) using the far RHS form of Eq. 4 and expressions for each half of the 
hemisphere located above the target surface. This approximation assumes that the variation in the azimuthal 
direction varies smoothly between values measured in the plane formed by the target normal and the ion beam. 
If the differential sputter yield data are curve-fit using Eqs. (6) and (7) for the front and back halves of the 
hemisphere located above the target., respectively, the total yield can be found analytically using Eq. (8). 

y(a) = A 6 cos 6 (a) + A 5 cos 5 (a) + A 4 cos 4 (a) + A 3 cos 3 (a) + A 2 cos 2 (a) + A 1 cos(a) ( 6 ) 

y(a) = B 6 cos 6 (a) + B 5 cos 5 (a)+ B 4 cos 4 (a) + B 3 cos 3 (a) + B 2 cos 2 (a) + B 1 cos(a) (7) 


6 

Y=I 


_i n + 1 


,_in + r 


The polar plot of the differential sputter yield data for normal (0°) and 30° incident xenon ions shown in 
Fig. 2c reveals dramatic differences between the two data sets more readily than those shown in Figs. 2a and 2b. 
The lighter gray lines shown in Fig. 2a and 2b that are labeled “Cosine Distribution” correspond to perfectly 
diffuse emission of atoms from the Ti target normalized to the total yield values listed in the figure legends. The 
cosine distributions were calculated using values of total yield found from Eq. (8) and the diffuse emission 
equation presented by Mahan: 11 

/ x Y cos(a) 

y(a) = — . (9) 

71 

Differential sputter yield data that lie below the cosine distribution at small polar angles and above at high polar 
angles are said to display under-cosine behavior (see Fig. 2a for example). The opposite condition of high 
differential sputter yields at small polar angles is referred to as over-cosine behavior (see the forward distribution 
shown in Fig. 2b for example). These same conditions are shown in Fig. 2c as a flattened circular shape for the 
1 keV, normal incidence ion beam (under-cosine) and a forward-extended cigar shape for the 30° incident ion 
beam (over-cosine). It is noted that an over-cosine distribution can also be generated from normal-incidence ion 
beams when higher ion energies are used. The results section of this paper contains more detail on sputtering of 
titanium and also presents data obtained to date on molybdenum, Poco graphite, and flexible graphite (grafoil). 

EXPERIMENTAL RESULTS 


In addition to the 1 keV data presented in Fig. 2, we have measured differential sputter yields of Ti at 
xenon ion energies of 0.75, 5, and 10 keV, and Fig. 3 contains plots of these data measured at normal incidence. 
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Figure 4 contains similar plots, but at an ion incidence angle of 30°. Several conclusions can be drawn from the 
data collected to date on Ti including (1) normal incident xenon ions at 1 keV and less sputter Ti into an under- 
cosine distribution, (2) higher ion energies of 5 and 10 keV at normal incidence cause over-cosine sputtering 
behavior, (3) significant forward sputtering occurs under oblique ion bombardment at low ion energies (0.75 and 
1 keV), but higher ion energies (-10 keV) result in more cosine-like behavior. 

Table 1 contains a list of the total yield calculated from curve fits to data like those shown in Figs. 2, 3, 
and 4. The front (A n ) and back (B n ) coefficients are also provided in Table 1, and they can be used in Eqs. (6) 
and (7) to calculate differential sputter yield distributions in the front and back halves of the hemisphere above 
the target (see Fig. 1). It is noted that all normal incidence data have been assumed to be symmetric and were fit 
in a way that forced the A n and B n coefficients to be equal to one another. A comparison to data collected by 
Rosenberg and Wehner 12 and a curve-fit suggested by Yamamura and Tawara 13 is shown in Fig. 5a. The data 
we collected fall -20% below the curve in Fig. 5a at low ion energy (0.75 and 1 keV) and -20% above the curve 
at high energies (5 and 10 keV). In addition to data on Ti, Table 1 also contains data collected to date on 
molybdenum, Poco graphite, and flexible graphite (grafoil). Figure 5b shows comparisons between our data on 
Mo at normal incidence and values obtained by Rosenberg and Wehner. 12 Our data were observed to fall 
slightly below the curve fit suggested by Yamamura and Tawara, 13 which was derived from data collected by 
several different researchers in addition to Rosenberg and Wehner. 

Differential sputter yield data are presented in Fig. 6 for xenon ion sputtering of molybdenum. The 
sample chosen for this test was a -15cm x 22 cm x 0.51 mm, arc-cast, cross-rolled molybdenum sheet that is 
typically used to fabricate grids for ion thrusters. The sample was tested in an as-received condition after light 
cleaning with a commercial de-greasing formula and water rinse. Figure 6a shows polar plots of data collected 
at 1 keV of ion energy including the corresponding data point on Fig. 5b. The normal incidence data were 
observed to be non-symmetric with more material being sputtered into the left half of the hemisphere above the 
target. The sputtered distribution from the 30° incident ion beam was observed to forward sputter Mo, but in a 
two-lobe manner. Finally, the 70° incident ion beam was observed to generate a nearly symmetric cosine 
sputtered distribution, with a slight over-cosine characteristic in the back hemisphere direction and a slight bulge 
in the front hemisphere direction. 

Figure 6b contains polar plots corresponding to 0°, 30°, and 70° incident ions at 5 keV. The normal 
incidence data (along with the 30° and 70° data) were observed to be much more symmetric than the 1 keV 
distributions, and all of the sputtered distributions at 5 keV are observed to be over-cosine in shape. Similar 
results were observed for sputtering by 10 keV ions as shown in Fig. 6c, except for the 70° data, which show 
high cosine-like, back- sputtered concentrations. 

A limited amount of data were taken using an electronic-grade graphite plate (Poco Graphite AFX-5Q) 
as the sputter target, and a polar plot of the data collected at a xenon ion energy of 750 eV is shown in Fig. 7a for 
angles of incidence of 0°, 30°, and 70°. For this sample, the normal incident beam was observed to sputter 
material in a highly under-cosine manner. As indicated in Table 1, the sputter yield for normal incidence was 
calculated to be 0.29 atoms/ion, which compares well to a value of 0.27 measured by Deltschew et al. 14 using PT 
2120, an electronics grade graphite plate from Carbone Lorraine. Data for 30° incident ions show a substantial 
forward sputtered character, while the data for 70° incident ions display a less forward and more symmetric 
over-cosine distribution. Figure 7b contains polar plots of sputtered grafoil (flexible graphite sheet) that were 
measured at 1 keV of ion energy at incident angles of 0°, 15°, and 30°. In contrast to the graphite sheet data, all 
the grafoil data at 1 keV were observed to be under-cosine, and, in addition, calculated values of sputter yield 
shown in Table 1 varied very little with incident angle. Similar trends were observed for grafoil at an ion energy 
of 5 keV where the sputtered yield remained nearly constant as the incident angle was varied from 0° to 30°. 

The largest incident angle-induced variation in sputter yield for grafoil was observed at 10 keV. It is also noted 
that the normal incidence sputter yield for grafoil is not much less than Ti at 10 keV. In addition, considering 
possible contamination from binders used in the fabrication of grafoil, this material would not be a good choice 
for the beam dump in a test facility designed to evaluate high specific impulse ion thrusters. 
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CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


Xenon ion sputtering yields were measured over an energy range from 750 eV to 10 keV and an angle of 
incidence range of 0 to 70 for several materials used in ion thruster devices and test facilities. The sputtering 
yield measurement technique utilized a quartz crystal monitor that was swept at constant radial distance from a 
small target region where a high current density xenon ion beam was aimed. This technique appears to provide 
sputter yield data that are in good agreement with the available literature. A particularly useful feature of this 
technique is the ability to obtain differential sputter yield data quickly at relatively low cost using modest test 
facilities. It is argued that these data would be valuable to analysts who are attempting to model the wear profile 
of an ion optics system throughout a given mission. Recommendations for future work include (1) extending the 
measurements to lower ion energies and higher angles of incidence and (2) documenting the differential sputter 
yield over the entire hemisphere above the target. Other important future work will be to measure the sputtering 
behavior of advanced carbon-carbon composite materials that are being considered for ion optics systems in the 
near future. Finally, we believe that evaluation of existing (and development of new) analytical and numerical 
sputtering models that give differential sputter yields should be pursued in parallel with an experimental effort in 
an attempt to “fill in the gaps” that would occur in an experiment-only program. 
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Fig. 2a Differential sputter 
yield measurements for Xe 
ions on a Ti target (normal 
incidence). 


Ion: 1 keV Xenon 
Incidence Angle: 0° 
Target: Titanium 
Yield: 0.57 atoms/ion 
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Fig. 2b Differential sputter 
yield measurements for Xe 
ions on a Ti target (30° 
incidence). 




Fig. 2c 


Polar plot of differential sputter yield measurements for Xe ions on a Ti Target (0° and 30° incidence). 
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Fig. 3a Differential sputter 
yield measurements for 750 eV 
Xe ions on a Ti target (normal 
incidence) 


Fig. 3b Differential sputter 
yield measurements for 5 keV 
Xe ions on a Ti target (normal 
incidence) 


Fig. 3c Differential sputter 
yield measurements for 1 0 keV 
Xe ions on a Ti target (normal 
incidence) 
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Fig. 4a Differential sputter 
yield measurements for 750 eV 
Xe ions on a Ti target (30° 
incidence) 



Fig. 4b Differential sputter 
yield measurements for 5 keV 
Xe ions on a Ti target (30° 
incidence) 



cos(a) 


Fig. 4c Differential sputter 
yield measurements for 1 0 keV 
Xe ions on a Ti target (30° 
incidence) 
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Fig. 5a Comparison of Ti sputter yield data to Rosenberg and Wehner 
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Fig. 5b Comparison of Mo sputter yield data to Rosenberg and Wehner 
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Fig. 6a Differential sputter 
yield measurements for 1 keV 
Xe ions on a Mo target ( for 0° 
30°, and 70° incidence) 
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Fig. 6c Differential sputter 
yield measurements for 1 0 keV 
Xe ions on a Mo target ( for 0°, 
30°, and 70° incidence) 
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Fig. 7a Differential sputter yield measurements for 750 eV Xe ions on a Poco graphite target 
( for 0°, 15°, and 30° incidence) 
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Fig. 7b Differential sputter yield measurements for 1 keV Xe ions on a flexible graphite (grafoil) target 
( for 0°, 15°, and 30° incidence) 
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